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ABSTRACT 

Two new transient testing techniques were evaluated; the centroid 
method developed by Kohlmayr and the time zero intercept technique. The 
zero intercept method was found to be the most promising of the two but 
is limited to values of N <C 2.5. The centroid technique can be used 
effectively when the value of N is less than 5.0. 

A heater system made of .001 inch diameter nichrome wire was designed 
and tested to determine its effect on the transient testing of matrix 
type heat exchangers. Because the design showed no improvement in the 
test results and was unreliable its use was discontinued. 
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English Letter Symbols 



NOMENCLATURE 



A Matrix tot*\l heat transfer area 

A Matrix uni.;n n free flow area 

c 

A r Matrix total ftor.tal area 

f r 

A Matrix solid cross-sectional area 

g 

available for thermal conduction 

a Fin thickness 

s 

b Flow passage perimeter (^A^) 

C f Fluid stream thermal capacity rate 

(mc f ) 

c^ Fluid specific heat 

C Matrix thermal capacity (W c ) 

S 8 S 

c Fluid specific heat at constant 

^ pressure 

c g Matrix material specific heat 

D Flow passage hydraulic diameter 

(4r h ) 

E Friction power per unit area 

G Flow stream mass velocity (m/A c ) 

g c Proportionality factor in Newton's 

Second Law 



g Normalized fluid temperature at 

Inlet to test section 

h Surface heat transfer coefficient 

for convection; heat transfer 
power per unit area per degree 
temperature difference 

1(g) Deviation from step 

k Fluid thermal conductivity 

^ Matrix thermal conductivity 



sq ft 
sq ft 
sq ft 
sq ft 

ft 

ft 

Btu/ (hr deg F) 

Btu/ (lbm deg F) 

Btu/deg F 
Btu/ (lbm deg F) 

Btu/(lbm deg F) 
ft 

hp/sq ft 

lbm/ (hr sq ft) 

32.2 (lbm ft)/(lbf 
sec^) 

dlmens ionless 
Btu/(hr sq ft deg F) 

dimensionless 
Btu/(hr sq ft deg F/ft) 
Btu/(hr sq ft deg F/ft) 
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L 


Total matrix flow length 


ft 


m 


Mass flow rate 


lbm/hr 


p 


Pressure 


Ibf/sq ft 


P 


Matrix porosity (A c /A ) 


dimensionless 


q 


Heat transfer rate 


Btu/hr 


R 


Gas constant (53.35-air) 


(ft lbf)/(lbm deg R) 


r h 


Hydraulic radius (A^L/A) 


ft 


t 


Temperature 


deg F 


t* 


Normalized temperature 


dimensionless 


u 


Flow velocity 


ft/sec 


V 

in 


Matrix volume 


cu ft 


w 

s 


Matrix mass 


Ibm 


w f 


Fluid mass in matrix 


Ibm 


X 


Distance along flow passage 
from the matrix inlet 


ft 


Z 


Reduced length (N x) 

L 


dimensionless 


Greek Letter 


Symbols 




ys 


Compactness (A/V ) 
m 


sq ft/cu ft 


p 


Ratio of orifice diameter to pipe 

diameter (d /d) 
o 


dimensionless 


A 


Difference or change (time, 
temperature, distance, etc.) 




e 


Time 


sec, hr 


A 


Free time 


dimensionless 


P 


Fluid viscosity 


lbm/hr ft 


M CENTO 


Centroid coordinate 


dimensionless 


P 


Density 


lbm/cu ft 
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Subscripts 



atm 


Local atmosphere 


ave 


Average 


f 


Fluid (gas, air) 


i 


Initial, inlet 


m 


Matrix, mean 


o 


At orifice 


s 


Solid (Matrix material), static 


STD 


Standard (temperature and pressure) 


X 


Local conditions 


1 


Inlet conditions (upstream of matrix and heaters) 


2 


Inlet conditions at matrix entrance 


3 


Exit conditions at matrix outlet 


Dimensionless 


Groupings 


f 


Fanning friction factor; ratio of wall shear 
stress to fluid dynamic head 


j 


2/3 

Colburn j-factor (N^Np^ ). This factor plotted vs 

Reynolds Number defines the surface heat transfer 
characteristics . 


A 


Longitudinal heat conduction parameter for solid 

material (k A /mLc-) 
s s f 


r 


Time parameter (hAG/W^c^) 
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1 . 



Introduction. 



The transient test facility at the Naval Postgraduate School (NPS) 
has been in operation for several years. Two techniques have been used 
at the NPS facility for the determination of heat transfer data. One of 
these is the maximum slope technique developed by Locke [13] which uti- 
lizes the temperature-time response curve of the fluid leaving the matrix 
after a step change in the fluid inlet temperature. Howard [6] extended 
this technique to include the effects of longitudinal conduction. The 
other technique is the cyclic technique developed by Bell and Katz [4]. 

In both of these techniques the heat transfer parameter being sought is 
N^^, which is a dimensionless heat transfer parameter equal to the ratio 
of the convective heat transfer rate from a solid to the heating capacity 
rate of an adjacent fluid, i.e., 

tia - 

where : 




= unit conductance for convection heat 
transfer (BTU/hr sq ft deg F) 

=» total heat transfer area of solid (sq ft) 

- mass flow rate of fluid (lbm/hr) 

* specific heat of fluid (Btu/lbm deg F) 

It is known that the maximum slope technique is unreliable at values 

of N less than 3,5, due to the large errors in N associated with 
tu * 6 tu 

errors in the determination of maximum slope. Furthermore, it has been 
noted experimentally that for values of N^ near 2.0, the temperature 
response curve of a fluid displayed its maximum slope at approximately 
time zero on the trace, and that this value of maximum slope was not 
suited for determining N 



A 

A 

rh 
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It was not until the work done by Kohlmayr [9, 10], which gave the 
exact analytical solution to the single-blow problem based on Hausen* s 
[5] mathematical model, that the reason for these errors was fully under- 
stood. Later work by Kohlmayr [11, 12] demonstrated ways in which to 
handle this problem as well as the development of a new technique which 
might be used in place of the maximum slope technique for values of N 
<5. Kohlmayr also developed a means to handle the single-blow problem 
when other than a step change in the inlet fluid temperature was made. 

The purpose of this thesis was first to determine experimentally 
the actual inlet temperature response of the NPS facility and to use this 
known temperature response in the manner suggested by Kohlmayr. Addi- 
tionally, a new heater system was made to try to more closely approxi- 
mate a step change in the fluid inlet temperature so that the physical 
process and mathematical model might more closely resemble one another, 
thereby improving the results using the maximum slope technique. 

Another technique, the "time zero intercept technique,'* used by 
Wheeler [17] was investigated. 
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2 . 



Summary of Theory, 



A, Background 

The single-blow transient technique which is used to determine 
heat transfer data for a porous solid originally used for its mathe- 
matical model Hausen } s [5] partial differential equation system. The 
method involved comparing the recorded exit temperature of a fluid passing 
through a porous solid which had previously undergone a step change in its 
temperature, with a computed response curve based on the solution to Hausen^ 
equations . 

Locke [13] has shown that there exists a unique relationship 

between the maximum slope of the response curves and the number of heat 

transfer units, N . However, Hausen f s model did not include the effects 
tu 

of longitudinal conduction, Howard [6], by the use of a finite differ- 
ence technique with the digital computer, included the effects of longi- 
tudinal conduction. 

It should be noted here that in all of Kohlmayr^ work, which 
will be discussed later, the effects of longitudinal conduction are not 
considered . 

B. Theory 

The basic assumptions in the single-blow problem are: 

(1) Properties of the fluid are temperature independent 

(2) Fluid flow is steady 

(3) The porous solid is homogeneous 

(4) The thermal conductivity of both solid and fluid is 
infinite in the direction perpendicular to flow 

(5) Thermal conductivity of the solid is zero in the 
direction of flow 
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Hausen's original differential equation system [5] was based on 
the energy balance between a fluid passing through a porous solid and the 



solid, The system of differential equations which resulted from this 
development are as follows: 

_ ts 

dj£_s ^3 = 

2r 



in which a "reduced length" variable, 

=■ - -/i/hc - // * 

2 ^ l 

and a dimensionless time parameter, 






M. 

l/\/c C<- 



e - 






% 



*S ^5 
have been used. 



m 



Ia // 






In the above equations: 

~ts = temperature of the solid (deg F) 

= temperature of the fluid (deg F) 

0 * time (hr) 

= specific heat of the solid (Btu/lbm deg F) 

Wi = mass of the solid (lbm) 

OC = distance along flow passage measured from inlet (ft) 

L ~ total length of solid (ft) 

K = mass of fluid entrained in solid (lbm) 

Kohlmayr [11] has modified Hausen's original equations by introducing a 
new dimensionless time variable called "free time," 



ll = mJk. g 

W'sX's 



Ws sC s L- 



(1) 
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X 



/VS 



'tvL 



( 2 ) 



For convenience Kohlmayr introduced two constants. 



= 02£c 



K A 



and 



z 5 



lA/ s ^< 



which are fixed for any particular experiment. For most practical pur- 
poses,^<f5 ^0 and ^0^0 . 

Hausen 8 s modified equations are then given as 

g)^/[ £ >/A + £ 5 ( 5 ,/^) 

(z,A) x t%(z.A) 

where t* and t* are the temperatures of the solid and fluid respec- 
tively* which have arbitrarily been normalized. The above equations are 
subject to the boundary and initial conditions: 

at H - 0 tf &/*-) = cZZ/V < 5 > 

and 

at/^-0 t* (Z,JLL) = / (6) 






- (it, A<-) 



(3) 



(4) 



where g ( yLL ) represents the normalized, time-dependent fluid tempera- 
ture at the inlet cross-section. Kohlmayr solved equations (3) through 
(6) by means of a double Laplace Transform [10]. 

The results of Kohlmayr 8 s solution are: 

iSM = e~ z [^h e^ (A ~ y) ^v)oiV\ 



+ 1-6 6 1 ^ cLy\ 



(7) 
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is an entire family of functions of the order K 0 J 



1 9 e a 



CO 



ZXX) = 




X 



MCn-wO\ 



n = o 



Evaluated at 

z.\ = N L 

IX-L 

equation (7) becomes 



/4. 



+ 



- /-<? [/- ^ 
A: s, ^ ' f/ - 



( 8 ) 



(9) 



By differentiation of equation (9) with respect to yLL , 

to(W) = eM&P + Kl z/W 









(10) 



A 

+ /C 



C -H, (^>-v;)l e <A 



With these results Kohlmayr [12] then showed that for other 
than a step change in inlet temperature, both the maximum slope and M free 
time" at a given value of N may be multi-valued. He further demonstrated 
that by knowing what the inlet response is, one can determine which of 
the values of maximum slope is valid and also how to determine what the 
error bounds are in the use of some known inlet temperature change. One 
of Kohlmayr ? s results is that the correct value of "relative" maximum 
slope occurring in the fluid exit temperature response will be that value 
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which occurred last in cases in which Kohlmayr^s technique is applicable. 

Furthermore, the maximum slope method is unstable for N <C 2.0, is 

tu 

singular at N = 2.0, and inaccurate for N between 2.0 and 3.0. [9,12] 

In view of the short-comings of the maximum slope technique, 
Kohlmayr developed an indirect curve-matching technique [11]. Some im- 
portant results of this development are summarized below. 

With the solution to the single-blow problem known, it is pos- 
sible to develop an indirect curve-matching technique based on the first 
moments of the fluid transient response curves. This method, known as 
the centroid method, involves the reducing of both the theoretical tempera- 
ture response anc * experimental response t *f ex p(^ tu > i nto 

two different single-valued functions based on the one parameter N . 

To generalize the problem further, a mapping functional was defined, 
based on the fluid inlet temperature change and the fluid's exit tempera- 
ture response such that 




where the following restrictions are imposed on ^4*; 

(1) y must be real, single-valued, continuous with 

respect to both ) anc * g(y^C ), and monotone with the parameter 

N • 
tu 

(2) 9^ (N ) must be monotone increasing with 1(g) = 

/oa ' tu 

J t the "deviation from step." 



(3) 






^ K 



, where K is some measure 



^4 A/tc^ 

for the maximum permissible amplification of errors. 

(4) For any N and any given deviation from step, ^4 (N^) 



should be insensitive with respect to local variations of g(y/X-) 
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(5) The evaluation of (N^) must be simple and 

straightforward. 

The functional chosen was the first moment of the difference 
between upstream and downstream fluid temperatures. 



oo 

f/* -y/Jfr- = 






■2 + M 



+ 






(U) 






Due to the difficulty involved with integrating equation (11) up to 
values of IM = Oo 9 a new functional was chosen and defined as 



*'4o r 










+ 



( 12 ) 



Some further restrictions imposed upon g(y^X) are: 

(1) g( ) niust be non-negative monotone decreasing 

(2) have initial value g( O ) =1 and 



(3) assume value zero for all free times which exceed 

the maximum permissible deviation from step: g( ) 

“ 0 for all n 21 I 

— max 

Therefore the moment functional was defined in terms of the centroid of 



the area under this difference curve: 






70 



A 






CENT£> 



— y o 



( 4 ^) - 4 ^ 



(13) 
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For the case of g {JUL ) = 0 (step change), the results of equation (13) 
are listed as Figure 1. 

In order to use the above results for any given inlet tempera- 
ture change g(^lL )» Kohlmayr arrived at the following empirical relation- 
ship: 



< 0) + • 99 1 ( 2 ] ± ■ 01 <“> 



This relationship was obtained by systematically applying the centroid 
method to many different upstream fluid temperature changes. 

Kohlmayr’s solution to the single-blow problem may also be used 
directly. Returning again to equations (9) and (1) 



= /- e v< ' 



1-fW 



+1 V 2 z, <? ' (i~y i-v))d~v 



where 



JM - m£r ° 



W/ -£/ 






W- Cc W s /c s i_ 

t^(N t ^,y^C ) can be evaluated at "free time," = 0, which would corres- 

pond to the time that an element of fluid which has undergone a change in 
temperature upstream arrives at the exit of the solid. Thus, 



' ‘ (f j, 

J t t+s 



-A/ 



i ~ Cj(Q) 



(15) 
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or 



I - %(%„<>) - [l-f(o)] 



(16) 



Taking the natural logarithm of both sides of the above equation, yields 



T tU 



= -M 






Jtn 



/- 



9 (c> 



(17) 



Therefore 



-JU 






H" Jin. 



-g(o) 



(18) 



Note that for a step change in the inlet temperature, i.e., g(ytC) = 0 
for yx> 0, that equation (18) reduces to. 



// = 

tu. 



/~t*lA/ ,0) 



(19) 



In order to utilize equation (19) directly it would be necessary to achieve 
a step change in the inlet fluid temperature, g(y£/-). This would occur 
at real time, Q =0. Then measure the temperature of that same element 
of fluid which had undergone the step change in temperature, as it left 
the porous solid. This would occur at "free time," LL. - 0. 
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Because of the finite response time of the different components 
of the experimental equipment* i.e., heaters, thermocouples and recorder, 
it is not possible to meet precisely either of the above requirements. 

In view of the above limitations it is necessary either to make some ap- 
proximations in the interpretation of the temperature response data, t*, 
or else devise a means by which the actual response data might be treated 
directly. Some of the difficulties involved with the second alternative 
are: Refer again to equation (9) 




/ - 



+ ■ J'kzstel/V-M cLiA 

o — j 



which may be stated as 



A/f.^ - |V ~ /('ll "f 



+ /'C 3, (A/ ( l cv-/o) e^: (y /U ~ \t - f(/i) ) ot v 

O 

and the expression for the inlet fluid temperature. 



where Q H is the equal to the heater time constant which was derived in 
Appendix C. 

One approach to using equation (9) would be to measure t* at 
some value of M free time'^ZX. > 0 such that g(yCC ) = 0 for which 

would be, 






A 1 , = “ O-i*/.) 



'6 

The difficulty here is 8 that in order to evaluate 

J M (t-ycv))dv 






23 



N must first be known* This would result in a direct curve matching 
tu 

technique* A second approach would be to measure t* at some time yl/^> 0 
such that /M- 

again N must first be known* In both the above methods if yOC is small 
the transient response of the sensing and recording equipment will still 
present problems. 

Because of the difficulties encountered with the above two 
methods, the first alternative, that of making some approximations in the 
interpretation of the measured response data, t*, will be considered. 

Since the response of the thermocouples is faster than that of the .003" 
diameter heaters, and approximately equal to that of the .001" heaters, 
see Appendix C, it might be reasonable to assume that the initial response 
recorded during a particular run is that response due to the heaters, see 
Figure 2B* Therefore, it is possible to extrapolate back to "free time" 

- 0, from the time on the response curve where the transient response 
of the heaters has died out and treat the intersection of the extrapolated 
curve and "free tim e,"^//4 = 0 as an actual step change in temperature. The 
above technique is referred to as the "zero intercept" technique. 
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3. Experimental Technique, 



Ihe existing apparatus at the NPS facility, see Figures 3 and 4, 
has been designed to conform to the idealizations required by Howard [5] 
for the use of his conduction parameter in the maximum slope technique* 
Howard's conduction parameter, ^ , is defined as 

\ As 

rh l jif 

where : 

As * thermal conductivity of the solid (Btu/hr sq ft deg F/ft) 

As * matrix solid cross-sectional area available for thermal 

conduction (sq ft) 

Howard's idealizations are as follows: 

(1) The fluid flow in the matrix is both steady and uniform 
in velocity and temperature at any cross section 

(2) The matrix thermal conductivity is finite in the direction 
parallel to fluid flow and infinite in the direction normal to flow 

(3) The matrix thermal conductivity is large in comparison to 
that of the contained fluid 

(4) The thermal properties of the fluid and matrix are constant 
and uniform 

(3) The convective heat transfer coefficient is some suitable 
average and remains constant 

(6) A step change in the temperature of the inlet fluid is im- 
posed at real time equal to zero. 

These idealizations result in making the fluid flow one dimensional. 
They in no way conflict with the restrictions imposed by Kohlmayr [11]. 

The requirement of uniform velocity and temperature profiles (1) 
is met by a specially designed entrance nozzle, flow straightening 
screens, and an even distribution of heater wires across the channel. 
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Piersall [15], using the equipment, verified that the velocity and tempera 
ture profiles were, indeed, uniform within reasonable limits. The im- 
portance of uniform velocity and temperature profiles was demonstrated 
by Wheeler [17] . 

The small temperature change across the matrix of about 20°F gives 
constant thermal properties of the fluid and matrix and a constant heat 
transfer coefficient as required by idealizations (5) and (6). 

The restriction that the temperature of the inlet fluid be subjected 
to a step change will be discussed later. 

The data required for the computation of the various heat transfer, 
fluid flow coefficients, and dimensionless parameters is as follows: 



p 

1 ATCA 


= atmospheric pressure (mm Hg) 


Po 


= orifice static pressure (inches 1^0) 


AS 


= pressure drop across orifice (inches 1^0) 


A $v 


= pressure drop across matrix (inches ^0) 


Ps 


= static inlet pressure at entrance to heaters 


to 


(inches ^0) 

= temperature of fluid at orifice (millivolts) 


cLo 


= orifice diameter 


p 

cs 


= ratio of orifice diameter to pipe diameter 
= chart speed (sec/inch) 


1 3 -~tz 


= downstream response 
= downstream - upstream response 



Pressures are measured with either a draft gage or a water manometer, 
depending upon the orifice-flow rate combination, except for atmospheric 
pressure which is measured by a mercury barometer. 

The temperature responses recorded are t^-t^, the difference between 
the inlet fluid temperature and the fluid temperature at the matrix exit 
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(used in the maximum slope and zero intercept techniques) and t^-t 
the difference between the fluid exit temperature and the matrix inlet 
fluid temperature (used in the centroid technique). For a more complete 
description of equipments, the reader is referred to Appendix A, See 
Figure 5 for the position of temperature and pressure measurements. 

A test run is accomplished by predetermining the necessary pressure 
drop across the orifice to achieve a desired flow rate and to determine 
the number of heaters necessary to achieve a 20°F temperature rise. Air 
is drawn through the apparatus and is controlled at the entrance to the 
turbocorapressor. When the desired flow rate is achieved the heaters are 
then energized and the heated air and test core are allowed to reach a 
steady state temperature, at which time pressure measurements are re- 
corded, Power to the heaters is then secured and recordings of t^-tj 
and t 3* t 2* both as functions of time are recorded on separate channels 
of the Brush recorder (see Figures 2A and B) . The temperature at the 
orifice (t^) is measured before and after a run to insure the same ambient 
air temperature for the run. 

After completion of the desired runs, the values needed to compute 
slope are taken from the recorded traces of t^-t^. The maximum value of 
slope is obtained visually with the aid of a straight edge. This informa- 
tion is included on the data sheet. The data sheet layout conforms to 
the data input section of a digital computer program [18], which reduces 
the data used in the maximum slope technique and calculates the parameter 
OC used in the centroid technique. For complete details on data re- 
duction the reader is referred to Appendix B at this time. 
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4. Description of Test Matrices. 

The two cores used in this experiment were both cores which had 
been tested previously and for which there was a good deal of data avail- 
able. 

One matrix, the Cercor T20-38, previously tested by Howard [7], is 
a ceramic type-core having a low thermal conductivity, and was used to 
eliminate any errors that might be introduced by longitudinal conduction 
which might have invalidated results using the centroid and zero inter- 
cept techniques. The other core tested, which is a stainless steel plate- 
fin type matrix. Solar 4, had been tested several times [15], [16], [18] 
and was used so that comparisons could be made when longitudinal conduc- 
tion was a factor. 

Further information on core geometries and properties is shown in 
Figures 6 and 7. 
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5, Presentation of Results, 

For each matrix tested, the heat transfer and flow friction character- 
istics have been computed*, The computed results are shown tabulated in 

Tables !-!!• The values of N have been plotted for the three methods 

tu 

used in this experiment,, Figures 8A and 9 compare values of N based on 
the maximum slope and centroidal technique for cores T20-38 and Solar 4. 
Figures 8B and 10 compare values of N based on the maximum slope and zero 
intercept technique for the same two cores*, Figure 11 is a comparison of 
N values computed by the maximum slope technique for the two heater 
systems. Core Solar 4 was the only core tested with the new .001 inch 
diameter wire heaters. 

The heater response curves are shown in Figure 12. See Appendix C 
for complete details used in determining these curves. 
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6 . 



Discussion of Results 



In order to insure that the technique used by the author in this 
experiment was correct , the results of the single-blow test of this experi- 
ment were compared with the results of previous tests of the same cores. 

Considerable difficulty was experienced in the use of the new (.001 
inch diameter) heater system. The heaters burned out frequently, thereby 
limiting the range of Reynolds Numbers over which the core was tested and 
also the number of cores actually tested. The exact reason for the burn- 
out is not known. It was initially thought to be due to surges in cur- 
rent caused by the use of a resistor bank in the first experiments. How- 
ever, when the resistor bank was removed and voltage to the heaters control- 
led by the motor generator set rheostat, the system failed again. In all, 
there were three different failures, all occurring at different flow rates 
and with a different number of heaters in use at the time of failure. 

An examination of Figure 11 shows that in the high Reynolds number 

ranges there is little difference in the value of N based on the maximum 

tu 

slope data for the two different heater systems. Furthermore, this is in 

a range where the errors in N due to errors in the maximum slope techni- 

tu 

que are greatest. Without actually applying Kohlmayr’s equations in the 
"extension of the maximum slope" technique [17], it is impossible to get 
a quantitative value for the error in N due to deviations from the step 
change. However, by the use of his curves based on the "deviation from 
step" (previously defined as 1(g) and equal to the area under the inlet- 
temperature response curve) which are presented here as Figure 13, it can 
be seen that for I < . 100 the error in slope. Am which is the difference 
in the maximum slope due to a step change in the inlet fluid minus the 
maximum slope due to a non-step change in the inlet fluid, is less than 
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.01 It must be pointed out that Figure 13, is for an upstream tempera- 

2 

ture change which is quadratic, i.e., g (JUL ) = (1-yU/A) . The curves of 
Am versus for the NFS facility would be slightly different from 

those presented in Figure 13 V for the : s associated with the NPS 

facility*s heater systems are exponential. However, for low values of 
I (I . 100) the error in maximum slope Am would be small regardless of 
the exact form of the Inlet temperature response. Since the maximum 1(g) 
resulting from either heater system is less than .03 it would not be 
expected for one to detect any difference of slope and consequently 
resulting from the use of either system. 

Again referring to Rohimayr f s curves, Figure 13, deviations from 
step changes in temperatures would cause large errors at lower flow rates, 
i.e , N^>5. However, since the 1(g) of the NPS facility decreases 
with decreasing flow rates (see Appendix C) it could be assumed that the 
errors in maximum slope due to 1(g) would be small. Furthermore, since 
the longitudinal conduction is greater at low flow rates its 1 effects 
cannot be neglected as they have in Kohlmayr's assumptions. 

In the investigation of the centroid and the zero intercept techni- 
ques the range of flow rates of both cores tested was the same and ranged 
from nr\ ~ 250 to DO vr 950 Ibm/hr. The reason for 950 lbm/hr being 

the maximum value tested was that at this flow rate all heaters were in 
use and in order to go to higher flow rates, an increase in voltage would 
have been necessary to achieve the same 20 deg F , temperature rise. This 
would ha-e changed rhe heater time constant, which is a function of both 
flow rate and the number of heaters in use. 

Investigation of the curves of Figure 8A and 9 reveals that for both 
cores tested there is very close agreement between the results predicted 
by the maximum slope technique, shown as the dotted line, and the results 



of the centroid technique for values of between 3*5 and 5* In fact 

there is even some overlapping of points in Figure 8A. For values of 

N below 3*5* the values of N based on the maximum slope technique 
tu tu 

begin to decrease from the predicted values whereas the values determined 
using the centroid technique followed closely the predicted values for 
the case of Solar 4 and increased slightly from the predicted values for 
core T20-38 * 

It can be seen from Figure 8B* which compares the values of N 
computed by the maximum slope and zero intercept method, for core T20-38, 
that the flow rates used were not high enough for the zero intercept 
method to be used to its best advantage* Only those runs made at flow 
rates greater than 600 lbm/hr had a clearly measurable zero intercept. 

Figure 10 again compared maximum slope N^ ff s with those evaluated 
by the zero intercept technique in this case for core Solar 4. This was 
done for both sets of heaters* It is important to note the close correla- 
tion between the values of N for the two heater systems. This will be 

tu 

covered more fully in Section 7* 
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7. Experimental Uncertainties* 

Various idealizations and boundary conditions have been imposed for 
the mathematical model of the physical experiment* Due to the fact that 
these idealizations and boundary conditions have not been precisely met, 
certain errors have been introduced. These errors are difficult to as- 
sign a numerical value to, and with the exception of the deviation from 
step temperature change previously mentioned and the effects of longi- 
tudinal conduction, will not be discussed* 

The experimental errors associated with the maximum slope technique 
used at the NFS facility have been covered in considerable detail, [2] [3] 
[15]; therefore only those errors introduced in this particular experi- 
ment will be discussed* 

In using Kohlmayrs centroidal technique there are several sources 
of possible error. First, Kohlmayr, in using Hausen p s mathematical model 
has neglected longitudinal conduction in the solid. Howard [6] has calcu- 
lated the errors in N associated with longitudinal conduction in using 
the maximum slope technique, but no such information is available for use 
with the centroidal technique, for the analytical solution to the single- 
blow problem, in this case, does not include the effects of longitudinal 
conduction. Therefore, one can only get a qualitative idea as to the ef- 
fects of these errors. 

From Howard c s curves of N versus maximum slope, Figures 14 and 15, 
for a given value of N the slope of a temperature response curve is 
different for different values of the conduction parameter A • There- 
fore, the response curves and consequently the centroid of the response 
curves would be different for different values of A . The difference in 

maximum slope due to A decreases as N decreases. For example at N = 

tu tu 

30 a difference in A of approximately .1 causes a 90% difference in 
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maximum slope*, whereas for N <T 5 the maximum error in maximum slope 

tu 

caused by a *1 difference in A is approximately 9%. It would seem 

reasonable to assume that provided Kohlmayr°s results were applied only 

in the low range (N^ <C 5) and longitudinal conduction was small (i.e. 

A . 1), conduction errors could be neglected. For this experiment 

the maximum value of A for either core at values of N <C 5 was less 

tu 

than ,022 . 



Another source of error is in the determination of the centroid of 
the area under the response curve. Two methods were used for this experi- 
ment , primarily as checks on each other. In both methods the critical 
point is the determination of where the cut-off is for a particular 
curve. The cut-off point is that point on the response curve which is one 
tenth of the maximum value. This point, defined as occurs at a 

point where the slope of the curve is slowly approaching zero, and an 
error of one millimeter in the ordinate can cause an error of several 
millimeters in the abscissa. For example, in the data taken with core 
Solar 4 when = 1.41, at (t^-tj) = 6mm, 0 * 13.0 sec., and for (t^-t^ 
=5, 0 =13.75 sec. This was for a typical trace in which = 

50.5mm. A sketch of the response curve (see Figure 16) will help to 
clarify the problem. 

The actual error in y^QENTD causec * ky an error in determining^^^ 
will depend on the value of N^. A typical example would be for 
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3.5. When computed by the manual technique (see Appendix B for details), 
it was found that a + 1 mm difference in ordinate (which is 2% and 
would be a maximum) caused an error of ^ 2% in It should be 

pointed out that these values are approximate, for there is some uncertain- 
ty associated with finding the intersection of the three centers of mass 
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lines, shown in Figure 16 and discussed in Appendix B. This uncertainty 
usually was less than + 27«. The same type of error associated with find- 
ing the cutoff point is also found in the computer technique of determining 
the centroid. In this method the computer solution required data points 
from the response curve at fixed intervals, which for this experiment were 
5ram intervals* Therefore, if the value of occurred at a point mid- 
way in any one 5ram increment, there could be as much as a 2% mm error in 
the cut-off. For an average 130 to 140 ram trace this is less than 2%. 

For comparative purposes, data runs were made using both techniques and 
the maximum deviation in centroid between the two methods was found to be 
4^, but the majority of deviations was less than 1%. 

Another difficulty discovered in analyzing the temperature trace was 
caused by the sensitivity of the thermocouples. During a particular run 
a change in the steady flow caused by a sudden draft in the laboratory 
caused a temperature deviation on the trace. This required some visual 
smoothing on the part of the author to extrapolate the actual response. 

This might cause errors when using an automated data reduction process, 
unless it included a good smoothing technique. 

An additional problem was encountered in the process of data reduc- 
tion. One of the assumptions made in using the centroid technique is that 
at time zero the temperature both upstream and downstream from the matrix 
is constant. It was found that this was not the case for this experiment, 
the reason being that there was some heating of the thermocouples located 
adjacent to the heaters due to radiation from the heaters. This meant 
that the zero reference point of the trace of (t^-tj) established before 
and after the run (heaters off) was not the same as that just prior to 
time zero when the heaters were deenergized (see Fig. 16). Therefore, 



35 



where the response curve should have resembled curve a of Fig. 16, it 
took the form of curve b. To compensate for this, point c was taken as 
time zero and curve a used to determine the centroid. It should be pointed 
out here that these figures are greatly exaggerated here for clarity and 
the error caused by this difference is felt to be negligible. Also, this 
difference was most evident at low flow rates, for the radiation effect 
was caused primarily by the heaters nearest the thermocouples. As the flow 
rates increased and more heaters were energized the percentage of heating 
due to radiation became less. 

So far no mention has been made of the error in N caused by an 
error Kohlmayr [11] has conducted a linear error analysis 

based on the approximation 

A ~ U- CENTO 

CENTO 



and introduced a relative error amplification factor 



K 



oL^t 



C.ENTP 



7^ 

c l/t-cEsrro 



or 

ANtm ~ «_ & M-CENTD 

y^CEAJTD 

With empirical results he has tabulated values of K versus an< * 

K versus N . 

tu 

As an example, using Kohlmayr°s curves when^^ E ^ D is measured too 
high by 2% and when K = -5.0, then N as obtained by this method would 
be 10% too low. 

Associated with the zero intercept technique are two distinct types 
of error. One results from the initial assumptions on which this technique 
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is based* These assumptions were discussed in Section 2. If errors are 



present in the results 9 they may be due to errors in the original assump- 
tions, It is not possible to get, at this time, a quantitative value for 
this type of error or in fact to even verify that the technique used here 
is valid. The only factor which would substantiate the technique °s vali- 
dity is that for both heater systems used in the experiment the results 
in terms of were nearly identical. The second type of error is that 
associated with the physical data reduction, primarily the correct extra- 
polation of the response curve back to time zero, see Figure 17. This 
may be seen in the following development: 




Mt^± AA^ = - Jyi 0~€* ± A if/) 





'-t't 

j± AM 



e 





or 




so that if 1-t* * .2 corresponding to a value of N 



tu 



1.609, and 




or 2% 



then 




or 1.25% 




or -2% 
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then 



0/98 




or -1.27o. 
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8 . 



Conclusions, 



1. The present design of the new ,001 inch diameter heaters was 
not found to be satisfactory due to their frequent and sometimes un- 
explained failures. No improvement of the results of the maximum slope 
and zero intercept tests was obtained using the smaller diameter wire 
heaters. No results were obtained with the new heaters for evaluation 
with the centroid technique 9 but it should be pointed out that for the 
centroid technique it is not necessary to have a step change in the in- 
let temperature so long as the actual change is known. 

2. It would be extremely difficult to try and compare the three 

techniques used in this experiment on the basis of test results ac- 
curacy. The best that one can do at present is to compare the various 
experimental results with the predicted values based on the maximum 
slope technique evaluated at high values of N (N > 5.0). In this 
respect the maximum slope technique is unreliable for N tu < ^ 3.5. The 
centroid technique appears to have good results for 5 to as low 

as .75 which was the minimum value tested in this experiment. For both 
the cores tested using the centroid technique the results were either the 
same as those predicted by the maximum slope technique or somewhat higher. 
The zero intercept method appears to be impractical for values of N^>2.5 
due to the larger differences associated with the logarithm of very small 
numbers. This might be better understood by looking at the governing 
equation in the zero intercept technique 

u. - — 

and noting that as the argument of the logarithmic term gets smaller the 
natural logarithm itself approaches negative infinity. Below the value of 
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N <T 2.5 the results were higher than predicted but were extremely con- 
tu 

sistent, that is, for different runs conducted at the same conditions 
the results were nearly identical* 

From the standpoint of ease of evaluating, the zero intercept techni- 
que is very fast and simple to apply. Furthermore, in the method used in 
this experiment the zero intercept technique appears to be insensitive 
to slight deviations from a step temperature change. On the other hand 
the results of the centroid technique rely entirely upon how well the so 
called "deviation from step" change of the fluid's upstream temperature 
is known. 

The centroid technique is much more time consuming to use and in 
view of the fact that its results are based on a considerable portion of 
the fluid temperature response trace it is more subject to errors caused 
by sudden fluctuations in the ambient conditions, for in the zero inter- 
cept technique only the very first part of the trace is of interest. 

Based upon the above conclusions the following table is recommended as 
a guide in determining which of the aforementioned techniques should be 
used for a given test range. 



N Range 

tu 


Technique 


3.5 < N 

tu 


Maximum Slope 


2.5 < N <3.5 
^ tu ^ 


Centroid Technique 


N <2.5 
tu 


Zero Intercept Technique 
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9 . 



Recommendations for Further Study, 



It is recommended that an attempt be made to determine both experi- 
mentally and analytically the effects of deviations from experimental 
assumptions such as constant fluid properties and the convective heat 
transfer coefficient on the results using the different techniques 
including the cyclic testing technique. 
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APPENDIX A 



Description of Equipment 

Air Supply, 

The working fluid for the experiment is air drawn through the test 
equipment by a 30 HP, multi-stage Spencer Turbo-Compressor, which is 
rated at 550 cfm operating on a 220 V a.c. power supply see Figure 3, 

Flow Measuring System* 

Flow measurement was accomplished with an ASME standard orifice 
section* Pressure taps were located d and d/2 diameters upstream and 
downstream respectively from the orifice. Thin concentric orifices with 
throat diameters of ,775, 1,232, 1,540, and 2,310 inches were used. 

Heater System, 

Two different heater systems were used. One system utilized a .0031 
inch diameter nichrome wire as the heating element, and the other .001 
inch nichrome wire. Both systems were designed to give the same nominal 
resistance for a given number of heater switches. 

The .0031 inch system, (Figure 18A) , consisted of 14 separate bakelite 
frames, each wound with two parallel-connected heater elements. Each 
pair of heaters was controlled by an individual switch, which in turn was 
wired in parallel with the other heater switches. All were controlled 
by one master switch. 

The .0010 inch system (Figure 18B) consisted of the same number of 
bakelite frames, with each frame containing six heater elements connected 
in parallel. The number of heaters and voltage could be varied to achieve 
approximately a 20®F. temperature change for any given flow rate. 

Matrix Holder and Test Section, 

The matrix holder and test section (Figure 4) were both constructed 
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of polyethylene plastic* The matrix holder is a drawer which slides into 
the test section, A removable frame on the side of the holder allows for 
the insertion of the matrix. The flow channel is 3-1/16 inches by 3-1/16 
inches, and can hold matrices up to 3 inches long. The matrix is held 
snugly in place by styrofoam insulation. On the downstream side of the 
holder is a plate containing the thermocouples used to measure t^. The 
test section into which the holder slides contains the heaters, pressure 
taps, and the thermocouple set downstream from the heaters, which measures 

V 

Inlet Cone and Flow Straightener 

This section was designed and tested by Piersall [15], and it pro- 
vided a uniform velocity profile to the air entering the matrix. 

Pressure Measuring System 

Pressure taps are located in the test section upstream and down- 
stream from the matrix holder, (see Figure 5). Two other taps are located 
at the orifice section. Each pressure tap is connected by flexible tub- 
ing to its corresponding manometer and/or draft gage. The following in- 
struments were used interchangeably, depending upon the flow rate: 

1. Ellison Draft Gage Company, 0-3 in. inclined gage 

2. Ellison Draft Gage Company, 10 in. manometer 

3. Ellison Draft Gage Company, 20 in. manometer 

4. Merriman Instrument Company, 120 in. manometer 

5. Precision Thermometer and Instrument Company, mercury 
barometer 

Temperature Measuring System 

Temperatures in the system are measured at four locations: (See 
Figure 5)inlet to the system ( t: ^ , downstream from the heaters (^) * 
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downstream from the matrix (t ) 5 and at the orifice (t ), 

3 0 

Temperature is measured by two different sets of 30 gage iron- 
constantan thermocouples 0 which were made by Traister [16] * Each set 
consists of 5 thermocouples connected in series* Each thermocouple is 
individually wrapped in teflon tape to prevent shorting, All ten are 
contained in an open-faced aluminum tube mounted in a frame at the exit 
of the inlet cone* The aluminum tube shielded against radiation from the 
heaters. 

Temperature t^ was measured by a set of five ,001 inch diameter 
iron-constantan thermocouples connected in series. The output of these 
was bucked against one of the sets measuring t^ so that the output of the 
two sets measured (t^-t^), 

t. 

Temperature t^ was measured in the same manner as t^ and the output 

bucked against the other set measuring t^ so that the output measured 

(t 0 -t-). The outputs (t 0 -t-) and (t^-t.) could be recorded separately or 
J 1 j 1 l 1 

bucked against each other to give (t^-t^)® The desired outputs were then 
led into an Astrodata Model 886 Wideband Differential D, C, Amplifier, 
where they were amplified 100s 1* From there the signal was led into a 
Brush Mark 280 Strip Chart Recorder, 

Temperature t^ was measured by a 30 gage copper-constantan thermo- 
couple referenced to an ice junction. The output was read on a Leeds & 
Northrup Millivolt Potentiometer, 

Heater Power, 

The power to the heaters was supplied from the 250 V, D, C, source 
in the laboratory. The actual voltage was controlled both at the supply 
panel and at the apparatus by means of a 50 SL wound resistor. 
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APPENDIX B 



Data Reduction Relationships 

(1) The following is a summary of the data reduction relationships 
used in calculating N by the maximum slope technique,, 

Geometry. 

Three geometric parameters may be used to define compact heat trans- 
fer surfaces. This allows for the comparison of different matrices. 

1. Hydraulic Diameter 

4x free flow area 



D._ = 4r, = , 

H h heat transfer area 



2. Porosity 



P m 



free flow area = c 



(B-l) 



(B-2) 



frontal area 



f r 



3. Area Compactness 



Q = heat transfer area = A 

matrix volume A, L 

fr 



(B-3) 



Dividing (B-2) by (B-3) 

r h = P/(S (b-a) 

Mass Rate of Fluid Flow 

The mass flow rate, m, is calculated from ASME Power Test Code [1] 
as modified by Murdock [14] by the following equation: 



m = 359 KaCh. yJWr 



(B-5) 



where 



c 



K = 



flow coefficient including velocity of 
approach 



V/ B * 

C = orofice coefficient of discharge [14] 

{3 = ratio of orifice diameter to pipe diameter 

d = orifice diameter in inches 
o 
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= thermal expansion factor 

$ =s P/RT the specific weight of fluid flowing^ assuming a 
perfect gas (lbf/ft3) 

Y - expansion factor 

P = absolute static pressure at orifice (Ibf/sq ft) 

R * gas constant for airs 53 0 35 (f t-lbf /lbm°R) 

T * absolute temperature at orifice (deg R) 

= pressure drop across the orifice in inches H^O 
Substituting the expressions for K and "^T in equation (B-5) yields: 

rh = 35-5 C- £ y [aEe 

V / -/3 V V fZ T 



From [1] Figo 40A 

Y = / ~ (o. V/ ~t O. 35 ft ¥ J A 



k = 1.4 for air 8 ratio 
Also, from [1], Fig. 38 




F « 1.0 
a 

P = (P -P /13. 6} (0.4912* 144) lbf/ft 2 
atm o 

P^tm ” local atmospheric pressure in inches Hg 

P = static pressure upstream of the orifice plate in inches Hg 
o 

Making the above substitution in (B-6) with the constants necessary to be 
dimensionally consistent yields: 



m = 58% 81 (C///-/5 1 ') cti 



!-(w 



■+ 




I 

J9.cz 




-to t-459 % 7 



(B-7) 
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Reynolds Number 



Reynolds Number is defined as: 

A/g e - A/ 

/X. 

where G. the mass flow velocity = m/A = m/pA,. 

c fr 

and JJ. is the fluid viscosity. 

Substituting: 

h/te = w &h B 



(B-8) 



(B-9) 



or 



y = 4 ‘mL / ftp. 



Maximum Slope 



(B-10) 



The maximum slope of the downstream cooling curve (t^-t^) is a unique 
function of N^ [13] and [6]. 












max 



^ , N^^ and ^ have been previously defined, but for convenience are 
restated : 



& } xM. and /) = 

M <s true, rintf L 

the new temperature introduced is t and is equal to the temperature 
of both the solid and the fluid prior to cooling. 

Therefore: 

^ /A^. c =■ (firiZf / W, 4 ) & 
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and 



oL Ct/A/, ) = (rn£t/ W $ yc s )de 



letting; 



and 



= mCj (Fluid stream thermal capacity rate) 



C “ W c (Matrix heat capacity) 
s as 

equation (B-ll) becomes*, 

<*(*/*£*.) = CCjCs) 

Furthermore 

^3 - tc_ _ 63 — £ f j 

*< ~ti t, -t; 



whose derivative iss 



cl (±?r^ 
\t, -±i , 



i cL It 3 -t , ) 



t 3-t 



combining equations (B-ll) and (B-12) yields 

^3 "’M 






-tt, 



oLfr/Ab) 



MAX 



_ Ka 






±,-t L 0 C 0 




(B-ll) 



(B-12) 



(B-13) 



I 
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From the above sketch of the downstream temperature response curve: 

= Y/X 
= d0 



do 

X/ chart speed 



MAX 



d(t 3 -t 1 ) 



Vh 



combining with the ratio of. 



matrix capacity 



Y 

A t 



flow stream capacity rate 



's (1/sec) 



and equation (B-13) 



c/i 



■Lt-Li 



d ( l K) 



f 'Y Chari SpeecL 

At x 



(B-14) 



_ fV\*X 



This value of maximum slope and /[ are then used to enter Table III 
or Figure 14 or 15 to get the corresponding value of N . 

(2) Centroidal Technique 

For the centroidal technique the recording trace of 
used. Two methods are actually involved in utilizing this trace. One 
method is to copy the trace physically with carbon paper on to a thin 
piece of cardboard. Then the trace copy is cut out, and by means of 
a r, plumb-line" attached to a pivot, which allows the trace to swing free, 
the line passing through the pivot point and the center of mass is determined. 
This is done for several points (at least 3) and the intersection of the 
different lines is the centroid of the area under the response curve. 
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These lines are shown in Figure 16. Ideally all lines will intersect at 
the same place * however * this Is not the case in actuality due to errors 
in manipulationo 

An alternate method is also used to find the centroid of the function 
This involves the use of a computer program* which when fed 
actual curve data (normalized) computes the centroid. This is done in 
the absence of an apparatus which could convert directly the temperature 
response to information useable by the digital computer* such as a paper 
tape puncher. 

Oncey^ ^ ^ is determined from the trace of necessar y 

to compute 1(g) in the manner described in Appendix Co This value of 1(g) 
is then used to determine^^^ • Recall Equation (14) 

JJ-CE^sro ~ + .99 !(<$■) ± .0/ 

or 

S. =Aw - - 99r ^ + - ol 

Once is determined* Kohlmayr 5 s curve of vs»JA 0 * Figure 1, 

may then be entered to give N 

tu 

(3) Zero Intercept Technique. 

To interpret the physical data* i.e.* the recorded trace of 
^ t 3* t l^ s> re P resents t^* in the manner described in Section 3, of 

this reports 

(i) Extrapolate by means of a French curve the response 
curve (t^-t^) from some position on the trace where the transient response 
of the heaters has died out* back to the vertical line passing through the 
point = 0. This point would be the position where the first change in 
temperature was detected by the thermocouples and is referred to as (t^-t^) 0 
or t f (N fcu ,0 )„ Figure 2B. 

(ii) Record the value of (t^-t^) at the intersection of 
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the time zero line and the extrapolated response curve, also record the 

value of (t^-t.) at the time just prior to the heaters being deenergerized. 

This temperature is referred to as (t 0 -t,) 

3 1 max 

(iii) N is then equal to, 

N tu = " ln(1 " (t 3' t 1^0 /(t 3' t l ) max ) 
which was developed in Section 3, as, 

4c = ■ t'-'t/WwO)) 
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APPENDIX C 



To determine the transient response of the NPS facility for differ- 
ent heater sets*, runs were conducted in which the temperature difference 
^2 ^1^ WaS recor< * e< * as a function of time. This temperature difference, 
which was previously called "’upstream temperature” is the difference in 
temperature of the air entering the test rig and the air leaving the 
heaters . 

Analytical determination of the heater time constant is straight- 
forward and presents no difficulty based strictly on the analysis of a 
cylinder in cross flow. The analytical determination of the thermocouple 
time constant is more difficult and requires some judgment. Since the 
thermocouples are constructed of .001 inch diameter wire, it is not pos- 
sible to describe precisely and mathematically the geometry of their 
junction. In the construction of the thermocouples the two different 
wires used crossed* at an angle approaching 180 degrees, and arc-welded 
at the junction. Therefore* the junction was neither spherical nor 
cylindrical (which is the case when the wires are butt-welded). For this 
reason the model of the junction was arbitrarily chosen to be spherical 
and calculations were based on a spherical junction .002 inches in dia- 
meter. On this basis it was determined that the time constant of the .003 
inch heaters was nearly three times greater than the time constant of the 
thermocouples* while the time constant of the .001 inch heaters was ap- 
proximately two- thirds that of the thermocouples. 

Since the response of the heater-thermocouple circuit would be the 
sum of at least two exponential terms* which would be difficult to 
separate experimentally* the terms have been treated here as one expon- 
ential term* which is a good approximation of the actual system response 
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if one neglects the initial points on the response curve. From the 
normalized temperature response curves the following empirical relations 
were developed? 







^2 ~t,) 






where 0 = 
6 h is 
in seconds 



time in seconds and 
the heater time constant 

(C-1) 




{^2 ”4 ) 
(tl-t,) ~M 






& locj & 



(C-2) 



6 




h\ 










(C-3) 



where the value of 



— 


Q 


log 


~ (t 2‘ t l ) 




^ t 2’ t l^ maX ' 



is the reciprocal of 



the slope of the temperature response curve. 

The heater time constant was then plotted against mass flow rate ( ) . 

From these curves the following relationships were obtained: 



log O/t = - a log rr\ + log c 

e„ - 



where a ^ c are (C-4) 

constants (C-5) 



Since the centroid technique makes use of the deviation from step defined 
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as 



Kf) = / d r ^ )d M- 



(C-6) 



'v 

and since ~ ® and 

then CJf^u) = fi(e) - e~ ° /& “ 



ct/j,' do 



1(9) 



-L e 



djjL — ^ 



(C-7) 



(C-8) 



where again 



- ihs£± 
Ka 



= m « constant (for any given 
core) 



Therefore to determine 1(g) for a particular run, enter the curve of 
vs. /'V? (Figure 12) with appropriate Yy\ . Note that this curve is 
independent of the matrix used in the run. Then multiply by the 

appropriate C< to obtain 1(g) used in centroid technique. Figure 12 
shows the experimental values of vs. />? plus the theoretical 

plots of heater and thermocouple transient response for both sets of 

0 

heaters . 
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Figure 2. Sample Recording of Temperature Responses 
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Figure 3. Experimental Setup 
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Figure 4, Matrix Holder and Test Section 
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Figure 5. Schematic of Test Apparatus 
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Figure 6, Geometric and Physical Properties of Core 

T 20-38 
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Figure 8. Core T20-38 Test Results 




Figure 9. 



Core Solar lj. Test Results 
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Figure 10. Core Solar !(. Test Results 
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Figure 11. Comparison of Maximum Slope Test Results of 
Core Solar 4> With .003" and .001" Heaters 
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"I jure 12. 



Heater Time Constant Versus Flowrate 
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Figure 13. Error in Maximum Slope Due to Deviation From Step 



NUMBER OF TRANSFER UNITS VERSUS MAXIMUM SLOPE 
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Figure 14. N^_ u Versus Maximum Slope With Conduction Parameter 

(Max. Slope, .3-2.05) 




Figure l£. N t Versus Maximum Slope With Conduction 
Parameter (Max. Slope . 3 - 1 » 7 ) 
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Finding U. From Temperature-Time Trace 
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Figure 17. Finding t„ From Temperature-Time Trace 
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Figure 18. Photo of .003 in. Wire and .001 in. Wire Heater Frames 
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*tu a3 a function of Maximum Slope and 
Longitudinal Conduction Parameter 
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N 2.5. The centroid technique can be used effectively when the value of N 

* tU i u c n tu 

is less than 5.0. 

A heater system made of .001 inch diameter nichrome wire was designed and 
tested to determine its effect on the transient testing of matrix type heat 
exchangers. Because the design showed no improvement in the test results 
and was unreliable its use was discontinued. 
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